The 3-D kinematics of water masers around the semiregular 

variable RT Virginis 
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ABSTRACT 

We report observations of water masers around the semiregular variable RT Virginis (RT 
Vir), which have been made with the Very Long Baseline Array (VLB A) of the National Radio 
Astronomy Observatory (NRAO) at five epochs, each separated by three weeks of time. We 
detected about 60 maser features at each epoch. Overall, 61 features, detected at least twice, 
were tracked by their radial velocities and proper motions. The 3-D maser kinematics exhibited 
a circumstellar envelope that is expanding roughly spherically with a velocity of ~8 km s . 
Asymmetries in both the spatial and velocity distributions of the maser features were found in 
the envelope, but less significant than that found in other semiregular variables. Systematic 
radial- velocity drifts of individual maser features were found with amplitudes of <2 km s yr . 
For one maser feature, we found a quadratic position shift with time along a straight line on 
the sky. This apparent motion indicates an acceleration with an amplitude of 33 km s yr , 
implying the passage of a shock wave driven by the stellar pulsation of RT Vir. The acceleration 
motion is likely seen only on the sky plane because of a large velocity gradient formed in the 
accelerating maser region. We estimated the distance to RT Vir to be about 220 pc on the basis 
of both the statistical parallax and model-fitting methods for the maser kinematics. 

Subject headings: masers — stars: individual (RT Vir) — stars:evolved, mass-loss, pulsation 
1. Introduction 
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Understanding the dynamics of mass loss flows 
in circumstellar envelopes of evolved stars is one of 
the most important areas of research in the over- 
all mass devolution of stars and the cycling of the 
interstellar medium. Fundamentally, material on 
the surface of an evolved star gets colder and forms 
dust while it is moving into interstellar space. The 
newly formed dust is accelerated by the stellar 
radiative pressure and forms an expanding enve- 
lope. However, such processes are complicated be- 
cause the related phenomena occur in a dynami- 
cally and physically unstable region. Maser emis- 
sion produced from silicon monoxide (SiO), water 
vapor (H2O), and hydroxyl (OH) molecules are 
commonly found in many circumstellar envelopes 
(e.g., Reid & Moran 1981; Elitzur 1992). Because 
of the compactness (down to 0.1 AU) of individual 
maser features, or probably physical gas clumps 
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(e.g., Colomer et al. 1992; Imai et al. 1997), they 
are good tracers for investigating the kinematics 
of the envelopes using very long baseline interfer- 
omctry (VLBI). Multi-epoch observations of SiO 
masers around Mira variables, which are located 
closer to a stellar surface than other masers, us- 
ing the Very Long Baseline Array (VLBA) suggest 
that circumstellar envelopes of Mira variables are 
not necessarily spherically symmetric, very likely 
due to anisotropic mass ejection on the stellar sur- 
face (e.g., Diamond & Kemball 1999). On the 
other hand, shock waves are also expected simul- 
taneously, which are driven by the periodic varia- 
tion in stellar radiative pressure to the dust in the 
envelopes. These shock waves caused by stellar 
pulsation are formed near the stellar surface and 
transported to larger distances from the star. 

The gas dynamics of circumstellar envelopes 
can also be studied with multi-epoch observa- 
tions of water masers with VLBI. The dynami- 
cal change of the envelope due to stellar pulsation 
is not so significant (^10%) in the near-vicinity 
of stellar surface, a so-called "radio atmosphere" 
and a "molecular atmosphere" (Reid & Mcnten 
1997). In fact, the amplitude of the expansion 
and contraction motions seen in SiO masers in 
this region is small (10-20%, Diamond & Kem- 
ball 1999). On the other hand, in the outer re- 
gion (T <1000 K) where water masers are ex- 
cited, dust formation starts and shock waves are 
enhanced due to dust-induced radiative pressure 
(e.g., Hofner, Feuchtinger, & Dorfi et al. 1995). 
A shock wave is expected to produce rapid veloc- 
ity changes (acceleration/deceleration) by up to 
10 km s _1 in front and in back of it, these will 
be directly detected as radial-velocity drifts and 
proper motions deviating from a constant velocity 
motion. The present VLBI technology enables us 
to detect such acceleration/deceleration motions 
with an accuracy of 10 microarcseconds (//as) in 
position and of 0.1 km s _1 in radial velocity. Such 
trials have recently started (e.g., Ishitsuka et al. 
2001, hereafter 101). Water masers are also good 
tools to directly estimate the distances to evolved 
stars using knowledge of their 3-D velocity vectors 
(radial velocities and proper motions), a method 
which does not rely on the standard distance lad- 
der methods (Marvel 1997). Distance estimation 
with water maser data will be applicable to dis- 
tance measurements for many evolved stars invis- 



ible at optical wavelengths after examining the re- 
liability and comparing with other methods. 

Here we present the 3-D motions and the radial- 
velocity drifts of water masers associated with the 
semiregular variable star RT Virginis (RT Vir), 
which have been measured from VLBA data. RT 
Vir is one of the brightest water maser sources 
(e.g., Bowers & Johnston 1994, hereafter BJ; Yates 
& Cohen 1994, hereafter YC; Richards et al. 
1999a, hereafter RCBY; Yates et al. 2000, here- 
after YRGB). The period of stellar pulsation of RT 
Vir has been estimated to be ~155 d with some 
irregularity (e.g., Kholopov et al. 1985; Etoka et 
al. 2001). Imai et al. (1997) (hereafter Paper I) 
found another pulsation period of 375 d using data 
obtained by the American Association of Variable 
Star Observation (AAVSO). Radial- velocity drifts 
of the water masers, with time, have been found 
from VLBI observations (Paper I), however, the 
origin of the drifts and the true maser kinemat- 
ics are not yet clear because of insufficient angular 
resolution and a small amount of proper motion 
data (c.f., RCBY; YRGB). 

Section 2 describes the VLBI observations and 
data reduction. Section 3 summarizes the revealed 
3-D kinematics of the water masers and the de- 
tection of acceleration of one maser feature in its 
proper motion. Section 4 discusses the dynamics 
of the circumstellar envelope of RT Vir, which ex- 
hibits acceleration. Measurements of the distance 
to RT Vir are also mentioned there. 

2. Observations and Data Reduction 

The monitoring observations of water masers 
around RT Vir with the VLBA have been made 
at five epochs during 1998 May-August, with 
a separation of 3 weeks between the successive 
two epochs. Table 1 summarizes the status of 
the observations. Each of the observations had 
a duration of 4 hrs including scans towards RT 
Vir and the calibrator 3C 273B for clock off- 
set and complex-bandpass calibration. The re- 
ceived signals were recorded with one base-band 
channel (BBC) with a center velocity at Vlsr = 
17.0 km s _1 and a bandwidth of 4 MHz in dual 
circular polarization mode. The correlated data 
were processed with 1024 velocity channels and a 
velocity spacing of 0.056 km s _1 in each channel 
at 22.24 GHz. 
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The procedures in data reduction using NRAO's 
AIPS and maser position measurement we have 
applied were almost the same as those presented 
in section 2.2 of Imai et al. (2000) and section 3.1 
of Imai, Deguchi, & Sasao (2002) (hereafter Paper 
III), respectively. In the present work, the com- 
plex bandpass characteristics were obtained from 
data on 3C 273B with an uncertainty of less than 
1° in phase. For amplitude calibration, we applied 
a template spectrum of the water maser emission 
obtained from the auto-correlation data. The ve- 
locity channel at Vlsr— 17.1 km s _1 was selected 
as reference for fringe fitting and self-calibration 
(see table 1). The naturally weighted visibility 
data created a synthesized beam of 0.41 mas x 
0.94 mas with a position angle of 7°. 6. The detec- 
tion limit was typically 200 mJy beam -1 at a 5-<r 
noise level in maps without bright maser emission. 
The positional accuracy for a maser feature was 
limited typically to 50 /zas mainly because of the 
extended structure of the feature. 

3. Results 

3.1. Proper Motions of Water Maser Fea- 
tures 

The water masers in RT Vir were spatially well 
resolved into individual maser features with the 
VLBA synthesized beam. Each maser feature is 
composed of a compact (<1 AU) bright part and 
an extended structure. As a result, about 60 maser 
features have been detected at every observation. 
Table 1 shows the numbers of maser features de- 
tected at individual epochs. Previous observations 
with larger synthesized beams (Paper I; RCBY; 
YRGB) may not be able to distinguish some of 
individual maser features because of the crowded 
distribution of maser features in the beam. 

We traced the position and radial velocity of 
each maser feature, which had been stable from 
one epoch to another within 1 km s _1 in radial ve- 
locity and 0.5 mas in position. We also checked the 
stability of spatial patterns of the maser feature 
(alignment of velocity components or spots, Paper 
III). A total of 60 maser features were detected at 
least twice and their relative proper motions and 
radial-velocity drifts were measured. Table 2 lists 
measured proper motions and radial- velocity drifts 
of maser features. The estimation of uncertainty 
in proper motion takes into account the apparent 



sizes of the maser features. Only 36 of the maser 
features were detected at least three times, indi- 
cating that individual maser features turn on and 
off on a time scale of typically 1-2 months. RCBY 
and YRGB reported maser proper motions of up 
to 3 mas in 70 days in the year 1996, the values 
were slightly larger than those measured in the 
present paper (typically 1-2 mas in 80 days in the 
year 1998). 

Figure 1 shows changes in positions and radial 
velocities of maser features detected five times. 
One can find jumps in positions and radial ve- 
locities within the above stability criteria in some 
maser features (e.g. RT Vir:I2002 20). These 
maser features were blended with nearby (closer 
than 0.5 mas) maser features, for which it was diffi- 
cult to exactly trace their velocities and positions. 
On the other hand, we found not only proper mo- 
tions that were well fit by a constant- velocity mo- 
tion but also exhibited linear radial-velocity drifts 
in some features. We also found an acceleration 
motion for a feature in its proper motion, which is 
described in detail in section 3.3. 

3.2. The 3-D Kinematics of the RT Vir 
flow 

Figure 2 shows the angular distribution and the 
3-D motions of maser features. The features ex- 
hibit an elongation in their spatial distribution 
and a radial- velocity gradient in the east-west di- 
rection, which have been visible in previous ob- 
servations (Bowers, Clausscn, & Johnston 1993; 
BJ; YC; Paper I; RCBY; YRGB). The 3-D maser 
kinematics clearly exhibits an expanding flow, not 
rotation. The kinematic flow looks clearer when 
including in the diagram the 60 maser features de- 
tected at least twice than when only including 36 
features detected at least three times. Because 
of the crowded distribution of maser features and 
very short lifetimes of maser features, misidenti- 
fication of feature motions should be taken into 
account especially for the features detected only 
twice. We made further analyses, however, using 
all of the proper motion data in the present paper. 

3.2.1. Analysis based on the velocity variance- 
covariance matrix 

We performed analysis based on diagonaliza- 
tion of the variance-covariance matrix (VVCM) 
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of our obtained maser velocity vectors (Bloemhof 
2000; 101). The velocity dispersions used in the 
WCM analysis are the quantities directly deter- 
mined by a proper motion measurement that lacks 
the absolute position reference. The WCM di- 
agonalization is done by obtaining eigenvectors 
and eigenvalues for the WCM, which correspond 
to the kinematic axes of the flow and velocity 
dispersions along the axes, respectively. Thus, 
the WCM analysis is an objective and model- 
independent method. The diagonalized WCM 
was as follows, 



0~XX 


0~yx 




0~xy 


a yy 








v. 



zy 



23.84 3.09 4.19 
3.09 24.11 7.17 
4.19 7.17 21.09 



15.03 
32.79 
21.22 



(1) 



where cry = Uji is a variance (i = j) or a co- 
variance (i 7^ j) of measured maser motions in 
the i- and j-axes (x, y, or z) in unit of km 2 s~ 2 . 
To estimate uncertainties of the obtained values, 
we also performed a Monte Carlo simulation gen- 
erating WCMs with artificial errors around the 
values obtained in the measurement. The eigen- 
vector corresponding to the largest eigenvalue (ve- 
locity dispersion) had an inclination of 24°.7±4°.6 
with respect to the sky plane and a position an- 
gle of 35°.8±14°.9, which is roughly parallel to 
both of the directions of the elongation and the 
velocity gradient mentioned above. This implies 
the bipolarity of the RT Vir flow, but is not so 
significant, since the velocity dispersion is roughly 
equal in all directions; a ratio of the eigenvalue 
was 2.47: 1.67: 1 (c.f., the ratio of 6.1: 2.0: 1 in R 
Crt, 101). 

3.2.2. Model fitting for the maser kinematics 

We also made a least-squares model-fitting 
analysis assuming a spherically expanding flow 
model, details of which have already been de- 
scribed in Imai et al. (2000) (c.f. Gwinn, Moran, 
& Reid 1992). We used weights proportional to 
the square of the accuracy of a measured proper 
motion. 

First, we adopted radial expansion motions of 
masers with independent speeds and estimated 



only a systemic bulk motion, or the motion of the 
star, and a location of the star as free parameters. 
Then, we also estimated a distance to RT Vir and 
the velocity field of the radial expansion. In the 
present paper, we assumed the speed of the radial 
expansion of a maser feature i, V exp (i), as a func- 
tion of the distance to a maser feature from the 
origin of the outflow, r,, which is expressed as 



W0 = V + V 1 



(2) 



where Vb is the intrinsic velocity at the stellar sur- 
face, V\ the velocity at a unit distance tq, a the 
power-law index indicating the apparent accelera- 
tion of the flow. We applied the fitting technique 
step-by-step, excluding maser features with unreli- 
ably large distances from the outflow origin (>0".3 
or 66 AU at 220 pc). Table 3 shows the best solu- 
tions obtained in the analyses. Figure 3 shows the 
distribution of features with expanding velocities 
obtained using equation (10) of Imai et al. (2000) 
versus distances from the estimated position of the 
central star. 

Through this procedure, we found that several 
maser features had negative expansion velocities 
or infall motions towards the star (see also figure 
2). Although some of the features were detected 
at least at three epochs, it is still unclear whether 
they were tracing actual physical motions. When 
including these infalling maser features, the best- 
fit function indicates a slow (V exp <5 km s _1 ) ex- 
pansion of the flow with marginal deceleration (fig- 
ure 3b). These are inconsistent with the sugges- 
tion that the RT Vir flow has an expansion velocity 
14 X p ~10 km s _1 in the water maser region and 
radial acceleration toward the outer OH maser re- 
gion (RCBY). The estimated position of the star 
is roughly at the center of the whole maser distri- 
bution but slightly biased toward clusters of the 
blue-shifted maser features located at the south- 
west of the whole region. The estimated position 
of the expansion center or the star is located close 
to the center of the "ring" found in the maser dis- 
tribution. No maser feature is located within 30 
mas ( 7 AU) from the star, indicating a "quenching 
zone" of water maser emission (YC). 

On the other hand, when excluding the maser 
features exhibiting infall, a higher expansion veloc- 



ity V e 



cxp 



;8 km s was obtained. The estimated 
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location of the star is very close to blue-shifted 
clusters of maser features. We could not find con- 
verging solutions in which the stellar position is 
close to the center of the ring mentioned above. 
Thus, the spatial distribution of the masers is ex- 
pected to be significantly asymmetric (see also fig- 
ure 2). This is consistent with the estimation by 
BJ. Maser features apparently spread out to 200 
mas (^45 AU) from the star (see also figure 3), 
which is larger by a factor of 2-3 than that pre- 
viously estimated (Bowers, Claussen, & Johnston 
1993; BJ; YC). The size of the maser- quenching 
zone is reduced to 10 mas (<~2 AU). The estimated 
systemic radial velocity of the star also has an off- 
set of ~3 km s _1 from the value adopted by above 
previous papers. 

Figure 4 shows the estimated 3-D spatio- 
kinematics of the water masers projected in three 
directions. Most of the maser features exhibit ra- 
dial expansion without rotation around the star at 
the diagram origin. Contrary to the asymmetric 
spatial distribution of maser features, asymmetry 
in the velocity field is not seen; maser features 
seem to have only a bias in the spatial distribu- 
tion. Thus, the kinematics of the RT Vir flow is 
well expressed by a radially expanding flow (c.f., 
Paper I; RCBY; YRGB) and the radial-velocity 
gradient in the east-west direction should be due 
to the bipolarity of the flow but with weak colli- 
mation. 

3.3. Acceleration motions found in water 
maser features 

We have measured radial- velocity drifts of indi- 
vidual maser features by measuring a radial veloc- 
ity at the brightness peak of each maser feature in 
the same manner as that in Paper III. Because we 
measured the velocities with accuracy better than 
a velocity channel spacing of 0.056 km s _1 , the 
uncertainties of the velocity drift rates were calcu- 
lated by assuming a measurement error to be equal 
to this value. The uncertainties were as small as 
0.3 km s _1 yr _1 in the best case. Table 2 gives 
the measured radial velocity drifts of maser fea- 
tures. Figure 1 presents time variations in radial 
velocities of individual maser features with longer 
lifetimes. 

The changes in radial velocity are usually less 
than the velocity widths of maser features (0.5- 
2.0 km s _1 , see table 2) and likely to be affected by 



temporal changes in spatial and velocity structures 
of maser features with time. Gwinn (1994) pointed 
out that the observed line width is likely affected 
by the subsonic bulk motions within the feature 
(<0.5 km s -1 ) and the different hypcrfinc transi- 
tions of the H 2 line (but negligible). Because the 
two effects will also affect the radial-velocity drift 
measurements, the velocity variations larger than 
0.5 km s -1 can be taken into account as possible 
accelerations. 

Nevertheless, some of the maser features changed 
their velocities with time at constant rates, or con- 
stant accelerations. Figure 5 shows the histogram 
of the measured radial-velocity rates. The drift 
rates of these features are equal to or smaller than 
1 km s _1 yr _1 and quite similar to those that have 
been commonly detected for 15 yrs with single- 
dish observations (Lekht et al. 1999). Maser 
features with larger drift rates (>2 km s _1 yr _1 ) 
exhibit jumps in velocities, which arc likely due 
to blending of a few very close by maser features 
as mentioned in section 3.1. Such large veloc- 
ity drifts were also found in some maser features 
detected in less than four epochs. Thus, radial- 
velocity drifts around 1 km s _1 yr~ 1 come from 
single maser features and may indicate accelera- 
tion/deceleration of the features. Previous VLBI 
observations (Paper I; 101) have also detected such 
systematic radial-velocity drifts in a small number 
of maser features. We did not find, however, a cor- 
relation between radial-velocity and acceleration 
(see figure 6). 

We have also looked for maser features ex- 
hibiting acceleration/deceleration motions in their 
proper motions. Similar to measurements of 
the radial-velocity drifts, we traced positions of 
brightest peaks of individual maser features. The 
method to define the feature position was de- 
scribed in Paper III. Note that flux-weighted po- 
sitions (e.g., Marvel 1997) are often affected by 
extended and weak velocity components. Figure 
7a shows the time variation in spatial structures of 
the position-reference feature (RT Vir 20021: 33) 
and the adjacent features. We find that the refer- 
ence feature has changed in spatial structure from 
one epoch to another. Most of the maser features 
had similar variation in their internal structures, 
which makes it more difficult to judge whether 
an apparent acceleration motion really reflects a 
true one or a physical motion of the maser cloud. 
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Because VLBI observations lose the absolute po- 
sition coordinate, we can not calibrate the change 
in spatial structure of the reference feature for the 
proper motion measurements. Assuming that this 
change is coming from the subsonic bulk motions 
within the feature mentioned above, an uncer- 
tainty of ~0.5 km s _1 (0.5 mas yr _1 ) likely con- 
taminates the measured relative proper motions 
of other maser features. 

Nevertheless, we found a maser feature candi- 
date that may exhibit a true acceleration motion, 
in which the feature moved by a distance much 
larger than the size of the feature itself. Figure 
7b shows the time variation of the candidate, RT 
Vir 20021: ^7, and the adjacent features. Figure 8 
shows the trajectory of the candidate and a model, 
which adopts a constant acceleration and is fitted 
well to the observed trajectory. The acceleration 
vector ([—18.5, 27.1] in units of km s _1 yr _1 ) has 
a position angle of (—34°, which is roughly paral- 
lel to the mean proper motion vector ([—2.2, 10.4] 
in unit of km s _1 , subtracting the mean proper 
motion shown in table 3) with a position angle of 
— 12° and the position vector ([—2, 46] in unit of 
mas) with a position angle of —2°. The measured 
acceleration motion vector is a relative one with 
respect to that of the reference feature, or we are 
measuring it in the frame in which the reference 
feature has a constant velocity motion. Expect- 
ing spherically symmetric acceleration motions in 
the features RT Vir 20021: 47 and 33 with respect 
to the central star, the relative acceleration vector 
would be closer to the proper motion vector of the 
feature RT Vir 20021: 47. Adopting this accelera- 
tion, the feature increased in speed from 1 km s _1 
to 16 km s _1 . 

Note that the apparent acceleration of the fea- 
ture was seen only on the sky plane (~33 km s _1 yr _1 ), 
but not along the line-of-sight (<0.1 km s _1 yr _1 ). 
Other maser features also exhibited apparently 
rapid acceleration motions on the sky plane but 
much smaller ones along the line-of-sight as men- 
tioned above (<3 km s _1 yr _1 ). 

4. Discussion 

4.1. Stellar position in the H2O maser dis- 
tribution 

The stellar position of RT Vir in the circmstel- 
lar envelope has been estimated in several manners 



based on the maser maps or position-velocity di- 
agrams. All of the results suggest that the star 
is close to the center of the H 2 maser distribu- 
tion or the middle between the blue-shifted and 
red-shifted clusters of maser features (e.g., Bow- 
ers, Claussen, & Johnston 1993; YRGB). YC and 
Bains et al. (2002) found a "quenching zone", 
within which no maser emission is found. Reid 
& Menten (1990) found that a radio continuum 
source from the semiregular variable W Hya is lo- 
cated at the center of a "ring" found in the H 2 
maser distribution. 

On the other hand, the dynamical center of the 
expanding flow estimated in the present paper has 
a large offset (~40 mas) over the uncertainty in 
the model fitting (see table 3). The stellar posi- 
tion can be examined by simultaneously making 
proper motion measurements such as those in the 
present paper and observations of the radio contin- 
uum source of RT Vir such as that made by Reid 
& Menten (1990). It is expected that the asymme- 
try of the 3-D distribution of maser features (see 
figure 4) derived a bias in the model fitting. Long- 
term observations, over several years, will reduce 
such distribution asymmetry. 

4.2. Distance to RT Vir 

The estimation of the distance to RT Vir has 
been made with several methods: a trigonomet- 
ric parallax with the HIPPARCOS satellite (140 
pc) and bolometric distance assuming a P-L re- 
lation (120-360 pc, e.q., BJ; YC; Yuasa & Unno 
1999). In the present paper, we made the distance 
estimation on the basis of the 3-D kinematics of 
the RT Vir water masers using both, the statis- 
tical parallax and the model-fitting methods. As 
described in section 3.2.1, the feasibility of these 
methods was examined by the objective analysis 
using VVCM, which shows no rotation or other 
peculiar motion in the velocity field (see also fig- 
ure 5). In the model-fitting method, the detail of 
which has already been described in section 3.2.2 
and table 3, we obtained distance values of 226±16 
pc in the analysis using only maser features ex- 
hibiting expansion. 

When including infalling features in the analy- 
sis, the distance value was 85±12 pc. Then we ob- 
tained an expansion velocity V cxp ^5 km s _1 (see 
figure 3), which is smaller than the value obtained 
by RCBY (V exp -10 km s" 1 ). The results of 
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RCBY suggested that the circumstellar envelope 
is accelerated gradually toward the outside, from 

2- 3 km s _1 in the SiO maser region closest to the 
star to 20-30 km s _1 in the OH maser region 
most distant from the star. It is expected that 
RCBY's results reflected a "steady" flow of RT 
Vir. Because, as expected in theoretical models 
(e.g. Hofner, Feuchtinger, & Dorfi et al. 1995), 
infall motions occur as a result of a time-variable 
phenomenon, or shock wave, the infalling features 
are needed to be excluded in the model assuming 
the steady flow. Thus, the above distance value is 
concluded to be underestimated. 

On the other hand, the statistical-parallax 
method (Genzel et al. 1981) gave distance val- 
ues of 237±29 pc (232±27 pc and 243±29 pc for 
the \x x -V z and the \x y -V z data, respectively) and 
214±28 pc (211±27 pc and 216±28 pc for the fi x - 
V z and liy-Vz data, respectively) when using all 

3- D motion data and only the data used in the 
former model fitting, respectively. The distance 
values obtained in the former model-fitting and 
the statistical parallax methods are quite con- 
sistent with each other, and the weighted-mean 
distance value of 220±30 pc is adopted for RT Vir 
in this paper. 

Note that the statistical parallax methods used 
only maser features with measured proper mo- 
tions. Adopting a radial velocity dispersion of 66 
maser features detected in the first epoch when 
the largest number of features were detected, a 
distance value of 279±32 pc was obtained. Thus, 
the radial velocity dispersion tends to have the 
larger value, with respect to the proper motion 
dispersion, when including maser features with- 
out measured proper motions. This implies that 
maser features without measured proper motions 
have the different kinematics from that of maser 
features with the proper motions. The distance 
estimation in this paper was the first successful 
one for a water maser source around a low-mass 
Mira-type star next to those around supergiants 
exhibiting many long-lived water maser features 
(Marvel 1997). 

4.3. Asymmetry of the mass-loss process 
in semiregular variables 

The morphology and the kinematics of circum- 
stellar envelopes (or mass-losing flows) of evolved 
stars, which are traced by the water maser kine- 



matics, sometimes exhibit clear bipolarity or sig- 
nificant asymmetry. Even if taking into account 
the limited physical conditions for maser excita- 
tion, such asymmetry is evident throughout the 
revealed 3-D maser kinematics. At the end of the 
AGB phase in stellar evolution, such asymmetry 
will be tightly related to formation of elongated 
planetary nebulae, some of which are created by 
stellar jets (e.g., Imai et al. 2002). The asym- 
metry, however, has been also found in semiregu- 
lar variables (e.g., Bowers, Claussen, & Johnston 
1993; BJ; YC; 101), which has been considered to 
be at the early- AGB phase. The present and pre- 
vious works (e.g., Marvel 1997: 101; Yates et al. 
2000) have proven that radial-velocity gradients 
seen in water maser distributions are due to the 
bipolarity, not the rotation of mass- losing flows. 
On the other hand, many Mira and semiregular 
variables, including RT Vir, do not clearly show 
the flow bipolarity yet, having major/minor axis 
ratios of <2. It is still difficult to discuss the re- 
lation among the asymmetry, mass-loss rate, ra- 
dius, and expansion velocity of a mass-losing flow 
due mainly to the limited number of sample stars 
available for data of the 3-D maser kinematics and 
for accurate distances. In this paper, we speculate 
about the weak asymmetry of the RT Vir flow. 

We re-estimated a mass-loss rate for RT Vir 
adopting the described spherically-expanding 
flow model. The density of material (hydro- 
gen molecules) in the flow is determined from 
predictions by a maser excitation theory as 10 6 
cm~ 3 < n H2 < 10 9 cm~ 3 (e.g., Cooke & Elitzur 
1985; Elitzur 1992). Adopting the radius of the 
flow in the water maser region, r ~44 AU, at 
d ~220 pc (distribution size of ~200 mas), an ex- 
pansion velocity, I4 xp ~8 km s _1 , we obtained a 
mass-loss rate of M ~ 1.5xl0 -7 (nn 2 /10 6 cm~ 3 ) 
Mq yr _1 , whose lower limit is roughly consistent 
with that estimated by BJ. Adopting the relation 
between a radius of the water maser distribution 
and a mass-loss rate found by Cooke & Elitzur 
(1985) and Lane et al. (1987), the mass-loss rate 
of RT Vir should be much larger than this value 
(c.f., 3xl0~ 6 M© yr-\ BJ). 

On the other hand, Mira-type stars with rela- 
tively large mass-loss rates and larger flow radii 
have symmetric kinematics (e.g., BJ). In OH/IR 
stars with larger mass-loss rates, 1612-MHz OH 
masers exhibit larger and symmetric distributions 
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(e.g., Chapman & Cohen 1986). It is considered 
that the intrinsic bipolarity or asymmetry of a 
mass-losing flow is obscured by the growing enve- 
lope as the mass-loss rate increases. In this case, 
RT Vir will be in the transition during an increase 
in its mass-loss rate. We have to take into ac- 
count, however, cases of supergiants, which have 
much higher mass-loss rates but most of which ex- 
hibit significant asymmetry. 

Note that the expansion velocity and the radius 
of a mass-losing flow determined seem to be time 
dependent. Water masers around R Aql are one 
of the examples (BJ). The bipolarity of the maser 
distribution was significant after the maximum of 
the light curve and vice versa. A similar varia- 
tion in the maser distribution has been confirmed 
for RT Vir, but when a stellar pulsation period 
of 375 d is adopted (Paper I). Adopting this ten- 
dency, water masers around RT Vir were observed 
around the light maximum, in which the maser 
distribution was larger than those observed with 
the VLA (Bowers, Claussen, & Johnston 1993; BJ) 
but smaller than those observed with the Japanese 
VLBI Network (J-Net, Paper I). If we could mea- 
sure maser proper motions with the J-Net data, 
the bipolarity would be exhibited more clearly. 
The time dependence of the bipolarity should be 
examined with monitoring observations for longer 
than a few cycles of the pulsation period. 

4.4. Origin of apparent acceleration mo- 
tions 

Long-term monitoring observations have showed 
that there are apparent radial-velocity drifts of 
water maser features that have occurred simul- 
taneously and exhibited the same drift direction 
(acceleration/deceleration) among several maser 
features (Lekht ct al. 1999). This implies a 
change of the velocity field in the mass-losing flow 
due to the stellar pulsation as suggested above. In 
the case of RT Vir, such acceleration/deceleration 
observed is likely excited by stellar pulsation with 
a period of ~375 d (Imai et al. 1997) rather than 
~155 d (e.g., Kholopov et al. 1985; Etoka et 
al. 2001). Such systematic velocity drifts may 
be detectable only with monitoring observations 
spanning at least one year. It is because there are 
small changes in the drift rates on this time scale 
and because it is sometimes difficult to remove 
data exhibiting velocity jumps due to blended fea- 



tures having intensity variation with time. In fact, 
in short duration (< a few months) the velocity 
drifts have not exhibited the same drift direction 
among maser features in both the previous and 
the present VLBI observations (see figure 6, also 
Paper I; 101). 

On the other hand, theoretical models adopt- 
ing shock waves driven by stellar pulsation predict 
that a shock wave creates acceleration/deceleration 
of a mass-losing flow and generates rapid velocity 
changes during the passage of the shock (a few 
months) (e.g., Hofncr, Fcuchtingcr, & Dorfi et al. 
1995). The velocity changes are expected to be 
^10 km s _1 and ^2 km s _1 at the inner (r ^5 
AU) and the outer (r ^50 AU) region of water 
maser excitation, respectively. The acceleration 
motion observed in a feature's proper motion, as 
mentioned in section 3.3, is quite consistent with 
the theoretical prediction. 

It has been a long-term issue whether change 
in the location of a maser feature is tracing actual 
physical movement of the gas in a clump rather 
than some kind of non-kinematic effect, such as 
traveling excitation phenomena or chance realign- 
ment of coherency paths through the masing gas. 
Recent observations, however, have strongly sug- 
gested the former in hydroxyl and water masers; 
all of the observations have shown that individ- 
ual maser features persist in both their spatial 
and radial- velocity patterns throughout their mo- 
tion for distances much larger than their own sizes 
(e.g., Blocmhof, Moran, & Reid 1996; Torrelles et 
al. 2001a, b; Paper III). The maser feature (RT 
Vir 12002:^7, see figure 7 and 8) exhibited a large 
acceleration on the sky, but simultaneously a large 
position drift of 2 mas, a stable radial velocity 
(variation less than 0.2 km s _1 ). This feature also 
exhibited a persistent V-shaped pattern formed by 
a cluster of maser spots, which have radial veloc- 
ities successively changing by 0.056 km s^ 1 from 
one spot to another. These strongly support that 
this maser feature traces the movement of a gas 
clump throughout five observation epochs. 

The larger acceleration on the sky plane 
(over 10 km s _1 ) than that in the line-of-sight 
(~1 km s _1 yr _1 ) is likely due to the beaming ef- 
fect of maser radiation. Taking into account the 
hydrodynamical treatment of a maser clump, an 
apparent acceleration motion of a brightness peak 

in the maser clump is expressed as, ~ V r • 
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where V r and r are the velocity field as a function 
of distance from the star, r, and the location in 
the clump, respectively. Assuming a velocity drift 
of 30 km s _1 yr _1 , a mean velocity of 8 km s _1 , a 
velocity gradient of ^£ ~4 k m g- 1 ma s _1 is ex- 
pected at 220 pc. This velocity gradient is larger 
than that observed in the above maser feature on 
the sky plane (~1 km s _1 mas -1 ) and that ex- 
pected from a steady accelerating envelope in the 
region of the H 2 masers estimated by Bains et 
al. (2002), 

dV AV 6 km s- 1 „ , , 
— — w — — ~ — — — — ~ 0.1 km s mas . 
dr Ar 20AU 

(3) 

If such a large velocity gradient is generated along 
the line-of-sight, the velocity-coherent path is too 
short to produce strong maser amplification occur- 
ring in a velocity width of less than 1- 2 km s . 
Adopting a velocity-coherent path 10 times as long 
as the observed feature size (~1 mas), an acceler- 
ation motion of less than 1-2 km s yr is ex- 
pected along the line-of-sight. Thus, only large 
acceleration along the sky plane is detectable. 

The present paper is the first approach to di- 
rectly detect the pulsation-driven shock waves, 
which should be identified in future works by find- 
ing simultaneously several maser features perform- 
ing such acceleration motions and by elucidating 
the relation between the occurrence of the accel- 
eration and the stellar light curve. 

5. Summary 

Our monitoring VLBA observations of water 
masers around RT Vir have revealed the 3-D kine- 
matics of water masers in more detail and de- 
tected the acceleration motions of water maser 
features. Because of the good angular resolu- 
tion (~1 mas) and time separations (~3 weeks in 
five epochs), we were able to measure 60 proper 
motions and radial-velocity drifts of maser fea- 
tures. Carefully dealing with maser features ex- 
hibiting infall towards the star, we obtained a 
spherically-expanding flow model from the 3-D 
maser kinematics. The estimated expansion ve- 
locity of ~8 km s _1 is consistent with that pre- 
viously suggested. The estimated outer radius of 
the maser distribution, ~45 AU, was 2-3 times 
as large as those previously estimated, except for 



those estimated in the J- Net observations. The 
estimated velocity field of the flow is roughly 
spherically symmetric, while the maser spatial 
and velocity distribution looks weakly asymmet- 
ric. Adopting the above flow model and taking 
into account a statistical parallax, we estimated 
the distance to RT Vir as to be 220±30 pc. We 
found radial-velocity drifts of maser features of 
<2 km s _1 and speculated that the drifts larger 
than >2 km s _1 are due to blending of a few maser 
features within a small region in space and veloc- 
ity. We also found an acceleration motion in the 
proper motion of a maser feature with a rate of 
33 km s _1 yr _1 and indicating radial acceleration. 
Such acceleration motions seen both in the line-of- 
sight and on the sky plane are explained by a sin- 
gle common phenomenon: passages of pulsation- 
driven shock waves in the circumstellar envelope. 
The possible variation in the maser expansion size 
and distribution asymmetry with time as well as 
observed acceleration motions should be examined 
in future observations with time intervals longer 
than a few cycles of the stellar pulsation. 
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Table 1: Status of the VLBA observations and position- reference feature in the proper-motion measurement 



Observation status Position-reference mascr feature (RT Vir: 12002 33) 











LSR 


Peak 


Position relative to the map origin 








Mascr 


Dopplcr 


intensity 








Date 




features 


velocity 


at Epoch 1 


x — Aq cos S 


y = AS 


Epoch 


(1998) 


UT span 


detected 


( km s" 1 ) 


(Jy beam -1 ) 


(mas) 


(mas) 


1 .... 


May 11 


03:00-07:00 


66 


17.14 


29.6 


83.885 


-6.922 


2 .... 


May 31 


01:30-05:30 


62 


17.16 


73.1 


0.000 


-0.069 


3 .... 


June 19 


00:00-04:00 


53 


17.10 


82.3 


-0.002 


-0.025 


4 .... 


July 12 


22:30-02:30 


50 


17.10 


84.2 


-0.015 


0.059 


5 .... 


August 1 


21:00-01:00 


53 


17.10 


63.5 


-0.034 


0.077 
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Table 2: Parameters of the water maser features identified by proper motion toward RT Vir 



Maser 
feature 
(RT Vir. 
12002) 


Offset' 5 '' 3 
(mas) 






Prop 
(m 


cr motion 11 
as yr" 1 ) 






Rad. motion 4 
(kms-l) 


RV dr 
(lm ,-1 


iff 1 ' 

yr - 




Peak intensity at five epochs 
(Jy beam- 1 ) 


116 

n a 


R.A. 


dccl. 


Mi 






(Tfly 


Vz 




Vz 




Ep. 1 


Ep. 2 


Ep. 3 


Ep 


. 1 


Ep. 5 


1 


















-2.98 






-7.57 






-0.02 






57.40 




110.33 






12.00 


Infall 


2 


90 


09 




53 


2 


02 


1 


02 


-0.26 


2 


18 


-7.53 






-0.82 


1 


38 








90 


08 


12.41 


Infall 


3 


63 


10 




82 




53 




60 


-10.58 




_ 4 


-7.17 




77 


1.74 




46 


24.65 




19.70 


48 


50 




Infall 


4 


83 


71 




99 


10 






25 


1.14 




98 


-6.89 




72 


-1.64 




31 


132.00 


47.40 


127.00 


81 


.51 


25.90 




5 


112 


87 


14 


68 


10 


95 




86 


-7.58 




59 


-6.39 




63 


0.24 




33 


4.09 


6.99 


4.56 






2.54 




6 


— 110 


.78 


26 


.95 


g 


.54 





29 


-2.08 





56 


-6.26 





39 


0.41 





.45 


14.20 


12.10 


8.32 


3 


30 




Infall 


7 


— 109 


.81 


2G 


.29 


_7 


.83 







-5.06 


1 


.35 


-6.18 





30 


-0.91 





74 


8.38 


5.35 


3.44 








Infall 


8 










1 






84 


-5.04 






-6.17 






0.73 






15.67 


12.46 










Infall 


9 


_ 112 


01 


9 9 


95 




62 


1 1 


24 


-0.31 




99 


-6.06 




51 


-2.39 


p 


45 


7.72 


5.66 


4.18 


2 


18 




Infall 


10 


115 


57 


16 


36 


"l4 


64 


1 


95 


-4.85 




77 


-5.99 




50 


0.87 


1 


38 








2 


24 


2.39 




11 


110 


97 


27 


40 


10 


35 




97 


-0.54 




59 


-5.84 




18 


-0.84 




74 


5.48 


6.22 


2.78 










12 


79 


21 


46 


60 




03 




33 


-10.06 




42 


-5.47 




86 


0.75 




45 


7.29 


9.04 


9.50 


12 


11 






13 


82 


47 


48 


46 




54 




24 


-9.20 




C3 


-5.30 




53 


0.82 


1 


45 


2.91 


1.98 










Infall 


14 


80 


38 


47 


83 




73 


2 


20 


-11.49 




29 


-5.22 




63 


-1.59 


1 


45 


4.31 


4.14 










Infall 


15 . . 


85 


89 


49 


'>(" 




63 




25 


-6.55 




05 


-5.20 




65 


0.82 




45 




7.86 


8.41 






4.98 




16 .... 


86 


21 


48 


79 


7 


20 




60 


-6.70 




24 


-5.14 




01 


0.53 


p 


31 


10.71 


14.30 


14.20 


8 


26 


10.70 




17 


— 87 


.32 


— 50 


13 


— 6 


.47 





.87 


-4.92 





79 


-4.57 





•1-1 


-0.02 


1 


38 








3 


66 


2.42 




IS 


— 87 


.59 


— 54 


.85 


— 6 


.53 





55 


-8.85 


1 


.26 


-4.03 





44 


-0.39 


Q 


4 6 


0.58 


0.51 







38 






19 


108 
















-5.74 






-4.03 






-0.86 










0.74 






2.45 




20 




35 


_ 5 g 


16 




^8 




38 


-0.13 




f 9 


-4.03 


II 


62 


1.78 


p 


31 


0.40 


0.76 


1.05 


1 


90 


0.21 




21 


105 


53 


49 






Or 


1 


82 


-8.48 






-3.97 




53 


-0.20 






0.59 


0.73 


0.68 










22 


87 


55 


53 


13 




2 3 




Vo 


-6.42 




45 


-3.97 




66 


0.22 




45 


1.42 


1.16 


0.76 





42 






23 


27 


73 


89 


82 


1 


73 






-7.42 




9 9 


-3.68 




31 


0.93 




45 




0.38 







38 


0.22 




24 


66 


42 


59 


16 


2 


38 


1 


52 


-10.97 




Of 


-3.66 
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0.33 


1 


45 


1.15 


0.72 












25 . . 


27 


■■") 
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1 


85 
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48 


-16.60 
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1.15 


1 


45 


0.42 
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26 


62 


84 


58 


57 


5 


38 





76 


-6.22 




58 
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26 


-0.06 


p 


31 
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64 
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27 
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1 Water maser features detected toward RT Vir. The feature is designated as RT Vir:I2002 N, where N is 
the ordinal source number given in this column (12002 stands for sources found by Imai et al. and listed in 
2002). 

2 Parameters at the first epoch while the feature was being detected. 

3 Relative value with respect to the location of the position-reference maser feature: RT Vir:I2002 34- 

4 Relative radial velocity with respect to the assumed systemic velocity of Vlsr= 18.2 km s _1 . 

5 Secular drift rate of a radial velocity. 

6 Infall: exhibiting a negative expansion velocity in the model fitting adopting only expanding features. 

7 Mean full velocity width of maser feature at half intensity. 

8 Uncertainty obtained when assuming a measurement error of a radial velocity, which equals to the channel 
spacing (0.056 km s _1 ). 
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Table 3: Best- fit models for the maser velocity field of the RT Vir outflow 



Including infalling features Only expanding features 



Parameter Step l 1 Step 2 2 Step l 1 Step 2 2 

Features finally used 46 48 

Offsets 

Velocity: 

V 0x 3 ( km s- 1 ) -1.9±0.7 -2.6±0.3 -3.2±1.4 -3.5±0.6 

V 0y 3 ( km s- 1 ) -0.7±0.5 -0.4±0.4 -4.3±1.0 -3.8±0.7 

Vq z 4 ( km s- 1 ) 5 0.7±0.3 5 2.5±0.5 

Position: 

x Q (mas) -39±6 -57±3 -11±9 -8±5 

y (mas) -13±5 -31±2 -14±6 -15±6 

Velocity field 

Radial outflow: 

Vq ( km s- 1 ) 6 -0.1±0.4 ... 6 6.6±0.6 

Vi ( km s-^rcsec"") ... 6 1.5±0.1 ... 6 2.0±1.5 

a 6 -0.40±0.03 ... 6 0.51±0.39 

Distance d (pc) 270 7 87±12 270 7 226±16 

RMS residual 7S* 3^02 6M 2^33 4.33 



1 Assuming independent expansion velocities of maser features. 

2 Asumming a common expansion velocity field expressed by equation (12) of Imai et al. (2000). 

3 Relative value with respect to the position-reference maser feature. 

4 Relative value with respect to VLsr= 18.2 km s _1 . 

5 Step 1 assumes the systemic radial velocity: Vq z = 0.0 km s _1 . 

6 The solution determines a radial outflow velocity Vexp(i) independently for each feature with a proper 
motion. 

7 Distance is completely covariant with the and Vexp(i) and cannot be determined: d =270 pc. 
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Fig. 1. — A sample of measured proper motions and the Doppler velocity drifts of maser features that had 
been detected in five epochs. The number added after "RT Vir:I2002" for each proper motion shows the 
assigned name. Solid lines in plots of proper motions show fit lines assuming constant velocity motions. The 
solid lines in the plots of the Doppler velocity drift show fit lines assuming constant acceleration motions. 
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Fig. 2. — The 3-D velocity field of water maser features around RT Vir. A 3-D velocity vector of a feature 
is indicated by a cone. The mean velocity vector of the maser features was subtracted from each of the 
observed velocity vectors, (a): the 3-D motions of 61 maser features, each of which had been detected at 
least twice, are shown. The plus symbols indicate the locations of the star, estimated by the model fitting 
using proper motions exhibiting both expansion and infall (pink) and only expansion (green), (b): Same as 
(a) but for 36 maser features, each of which had been detected at three epochs. 
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Fig. 3. — The expansion velocity of a maser feature against the distance from the position of the central star. 
These parameters were estimated after the model fitting. The expansion velocity was calculated by using 
the equation shown by Imai et al. (2000). A solid line shows the velocity field as a power-law function with 
respect to the distance from the star, (a): For the case when excluding infalling features, having negative 
expansion velocities, (b): For the case when including all maser features within a reasonable distance from 
the star (< 300 mas, corresponding to 66 AU at 220 pc). 
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Fig. 4. — Estimated 3-D positions and motions of water maser features in RT Vir. The positions and 
motions are with respect to those of the outflow's origin determined by a model fit. The position of the 
arrow indicates that of a maser feature. The direction and length of the arrow indicates the direction and 
the magnitude of the maser motion, respectively. Top: Front view (XY-plane) for the positions and motions 
plot. Bottom left: Top view (XZ-plane). Bottom right: East-side view (ZY-plane). 



1G 




Fig. 5. — Histogram of drift rates of maser features' radial velocities. Maser features are divided into two 
groups, features detected at all epochs and those detected at 3 or 4 epochs. Maser features detected only at 
two epochs have been excluded. 
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Fig. 6. — Dependence of the radial- velocity drift of a water maser feature on its location and radial velocity, 
(a): Against distance from the star (b) Against radial velocity with respect to the assumed stellar velocity 
(Vlsr=18-0 km s _1 ). No correlation is found between radial- velocity drifts and locations/radial velocities. 
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Fig. 7. — Spatial structures of maser features and 
their time variation in the two selected fields. The 
number added after "RT Vir:I2002" for each maser 
feature shows the assigned name. Each of the filled 
circles shows a velocity component (maser spot). 
Radial velocities of spots change by 0.056 km s _1 
from one spot to the adjacent spot, (a): Maser 
features around the map origin. The position- 
reference feature, RT Vir: 12002 33, is spatially 
fixed at the map origin to make measurements of 
maser proper motions, but will move from south 
to north so that it expands with respect to the 
star, (b): Maser features exhibiting acceleration 
motions. The maser feature on a black line, RT 
Vir:I2002 ^7, shows a clear constant acceleration 
motion and the development of its V-shaped struc- 
ture. The maser feature had kept its brightness 
peak close to the crease of the V-shaped structure. 
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Fig. 8. — The temporal position variation of the 
maser feature RT Vir:I2002 ^7, which apparently 
shows an acceleration motion. Solid lines indicate 
fits to the position variation in R.A. and decl. di- 
rections assuming a constant acceleration motion. 
The positions of the maser feature at individual 
epochs were defined to be at the brightness peak 
among the maser spots in the maser feature and 
measured with uncertainties less than 100 /ias. A 
vertical bar with the feature is an extension size 
of the feature, or a distribution size of the spots 
in the feature. 
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